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Abstract 

We develop a method to measure the CKM angle 7 without hadronic uncertainties from the 
analysis of — » K^D and K^D followed by singly Cabibbo-suppressed D decays to non CP- 
eigenstates, such as K^K*^ . This method utilizes the interference between b — > cus and b — > ucs 
decays, and we point out several attractive features of it. All the modes that need to be measured 
for this method are accessible in the present data. 
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Some of the theoretically cleanest determinations of the weak phase 7 rely on B — > DK 
and related decays 0, fl 3, Q (for definitions, see Q, Q). The original idea of Gronau 
and Wyler (GW) [J| was to measure two decay rates arising from b — > cus and b — > wcs 
amplitudes. By measuring the rate of a third decay that involves the interference between 
these two amplitudes, one can gain sensitivity to their relative phase, which is 7. Since all the 
quarks which appear in B — > DK decays have distinct flavors, the theoretical uncertainty 
arises only from higher order weak interaction effects (including, possibly, D — D mixing, 
which we discuss below). However, there are no penguin contributions to these decays. 

A practical difficulty of the GW method is that the amplitude ratio A(B~ — > D°K~)/ 
A(B~ — > D°K~) is expected to be small. As a result, the measurement of \A(B~ — > 
D°K~)\ using hadronic D decays is hampered by a significant contribution from the decay 
B~ — ► D°K~, followed by a doubly Cabibbo-suppressed decay of the D°. To avoid this 
problem, Atwood, Dunietz, and Soni (ADS) |s| proposed to study final states where Cabibbo- 
allowed and doubly Cabibbo-suppressed D decays interfere. Several other variants of the 
GW method have been proposed jf], 0| • An important point is that most of the methods 
require the measurements of very small rates, which have yet to be observed. (One exception 
is Ref. where relatively large rates are expected.) 

In this letter we propose to use singly Cabibbo-suppressed D decays to final states that 
are not CP eigenstates, and have sizable rates in both D° and D° decay. The use of such 
final states has been mentioned in Ref. j^; here we develop the details and point out the 
advantages. The simplest example is the final states K* + K~ and K*~K + . In the formalism 
below, we assume that all relevant decays are dominated by the standard model amplitudes. 1 
We define 

A B = A(B~ -> D°K-) , A B = A{B- -> D°K~) . (1) 



We shall further denote 



Ab 
A b 



r B e 



A, 



A } 



(2) 



where 5b is the relative strong phase between A B and A B , and we have neglected the 
deviation of the weak phase of Ab/Ab from —7, which is suppressed by four powers of the 
Cabibbo angle (A 4 ~ 2 x 10~ 3 ). Then the ratio of the CP conjugate decay amplitudes is 



1 New physics would have to compete with tree-level b — > ucs or c — > ssu decays, or give rise to D — D 
mixing near the present experimental limits to influence our determination of 7. 
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given by 



A(B + -> D°K^ 



A(B+ -> D°K + ) 
We denote the following .D decay amplitudes 
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and their ratio 
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Here we neglected the c —>■ u penguin contribution compared to the c — > ssu tree diagram, 
which is a very good approximation. Then the ratio of the CP conjugate amplitudes is 

A(D° — ► K*~K + ) „ 



r D e 



A(D° -> K*-K+ / 

With these definitions, the four amplitudes we are interested in are given by 
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K+(K*-K + ) D ] = \A B A D \ 
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K+(K* + K-) D ] = \A B A D \e 



rD + r B e 
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(7) 
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Of the unknowns in these equations, \Ad\ and have been measured in D decays and 
\A B \ was measured from the B~ — > D°K~ rate (and its conjugate) by reconstructing the D° 

nnn 

in flavor-specific decays [9|, [10|, ■ (While in practice, measuring \A B \ involves identifying 
the D° through its hadronic decay, which is not a pure flavor tag, this induces a negligible 
error.) For any given integrated luminosity in the future, the errors in the measurements of 
| Ad I) r Di and \A B \ will induce a smaller error in the measurement of 7 than the statistical 
error of measuring the decay rates corresponding to Eqs. (JJJ). 

This brings us to the key point: by measuring the rates of the four decays in Eqs. ( |7| ). 
one has four measurements for the remaining four unknowns: r B , 5 B , Sp, and 7. 

A simple analytic solution for 7 can only be obtained neglecting the terms proportional to 
r 2 B in the four decay rates corresponding to the amplitudes of Eqs. (J7J). In this approximation, 
we obtain 

2 (Rt + R 3 - 2) 2 - (R 2 + R 4 - 2r 2 D ) 2 



cos 7 



4 [(R, - l)(R 3 - 1) - (R 2 - r D ){R A - r D )\ 



(8) 
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where 

A[B-^K-(K'+K-) D ] 2 

Rl = AjC D ' (9) 

and similarly R2-4 are the squares of the "reduced amplitudes" corresponding to lines 2 — 4 

in Eqs. (|7|). Eq. (JHJ) illustrates that our method is sensitive to 7 even neglecting terms in 

the branching ratios proportional to r%. Although r B is expected to be small, r B should of 

course not be neglected when the experimental analysis is carried out numerically. On the 

other hand, r B is expected to be of order unity, and our method works best if rjj and 1 — r_n 

are both of order unity. This expectation is supported by the data, td = 0.73 ± 0.21 8| 

(the real uncertainty of rjj may already be smaller; here we assumed that the errors of the 

measured D° —>■ K* + K~ and D° — > K*~ K + rates are uncorrelated) . 

In principle, our method requires the analysis of only one type of final state with different 
charge assignments (such as K^K^). In practice, the sensitivity can be improved by 
considering several B and D decays of the type considered so far. When a different B decay 
mode is used, for example, B~ — > K*~ (K* + K~)d, four more measurements can be done, 
but only two new parameters are introduced, r' B and 5' B . (Here we assumed, as before, that 
\A' B \ of this B decay channel is measured.) When a different D decay mode is used, for 
example, B~ — > K~(p + tt~)d, only one new parameter is introduced, 5' D . (Here we assumed 
again that \A' D \ and r' D of this D decay channel are measured.) An especially interesting 
additional D decay mode is to CP eigenstates. In this case, two extra measurements are 
possible, but no new parameters are added, since r' D — 1 and 5' D = or n. 

Next we discuss how the sensitivity to 7 changes in some limiting cases. If ro = 1 then 
Eqs. (JIJ) become degenerate if either 5d = or 5b = 0, and 7 can no longer be extracted 
from this mode alone. If ro 7^ 1, then our method is sensitive to 7 independent of 5d and 
5b- However, if 5d = or 5b = then the sensitivity to 7 comes only from terms in the 
branching ratios proportional to r%. In this case Eqs. (J7J) are not degenerate, but both the 
numerator and the denominator of Eq. (JHJ), obtained neglecting r B terms, vanish due to 
I cos(5d + 5 B ) I = I cos(5d — 5 B )\. This indicates that the error in the determination of 7 
may become large if either of the strong phases is small. This potential difficulty may be 
eliminated using several decay modes, as long as there are two sizable and different strong 
phases. For example, even if 5b is small, the sensitivity of our analysis to 7 can still be large 
and not rely on the r B terms in the decay rates if we use two D decay modes with sizable 
strong phases. 



Throughout this analysis we assumed that |Ad| and rjj are known from D decays, but 
Sd is not. In the near future, it will be possible to measure 8d at a charm factory using CP- 
tagged D decays simplifying our 7 measurement. With some model dependence, So 
may also be measured at the B factories using a Dalitz plot analysis of the D decay (l5( . In 
this case, one typically assumes that the variation of 5r> over the Dalitz plot can be accounted 
for by using phases of Breit-Wigner resonances. With enough events, the validity of this 
assumption can be checked in the analysis. Both the charm and the B factory measurements 
of T£, and 8jj can be carried out as a function of the Dalitz plot variables. 

Measurements of 7 that depend only on one strong phase, S, are in general subject to an 
eight-fold discrete ambiguity, due to invariance of the observables under the three symmetry 
operations jl^ 

Sex ■ 7 -* 8 , 8 -> 7 , 

SWn : 7^-7, 8-* -8, (10) 

S-k ■ 7^7 + 7T, 5 ^ 5 + TT . 

In the modes we propose, the variation of 8d across the D decay Dalitz plot will be largely 
determined by the Breit-Wigner shape of the dominant resonance (such as the K*). As a 
result, the theoretical expressions for the decay rates are no longer invariant under S ex and 
Ssign, and the only ambiguity that is relevant for our method is the two-fold S n ambiguity Q. 

To compare the sensitivity to 7 of the different methods, we assume that the only small 
parameters are A and tb- The latter has been estimated assuming factorization as 

v ub v: 



^~A/2, (11) 



v a y;, 

where Nc = 3 is the number of colors. The accuracy of this estimate is expected to depend on 
the specific hadronic mode. Thus, for example, the numerical values of tb in B — > K(K*K)jj 
and in B — > K*(K*K)jj are expected to be different. Since the uncertainty of this estimate 
of Tb is large, it is not yet known whether is closer to A or to A 2 . The statistical 
significance of a CP asymmetry measurement scales roughly as the smallest amplitude that 
is needed in order to generate the asymmetry. Thus, to compare the methods we need 
to identify the smallest such amplitude in each of them. Compared to Ab, the smallest 
amplitude in the ADS method is of order min(A 2 , r#), while in the GW and in our methods 
it is of order r#A. This simple argument suggests that if r# < A then the ADS method 
has the largest sensitivity while if tb > A then it is one of the others. There are many 
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additional factors and experimental differences that will influence this comparison when 
the measurements are actually carried out. For example, since in our case measurements 
of doubly Cabibbo-suppressed D decays are not needed, the induced experimental error 
from D decay rates is expected to be the smallest. We conclude that the sensitivity of these 
methods are comparable and depend on yet unknown hadronic amplitudes and experimental 
details, and so all should be pursued. 

So far we have neglected D — D mixing in our analysis. Its effects on the GW and ADS 
methods were studied in Ref. [l^, and it is straightforward to generalize it to our case. We 
find that the leading effect on our determination of 7 is generically of order xd/^b and Ud/^b-, 
where Xd — Am D /r D and yo — Ar D /2T D . Since the present experimental bounds on xd 
and ud are at the few percent level, and even the standard model values may be not much 
smaller [16], D — D mixing gives rise to a theoretical error of order 10%. This is the largest 
theoretical uncertainty in our method at present, but it will be reduced as experiments yield 
tighter bounds on or measurements of xd and yn- (The leading sensitivity to D — D mixing 
in the GW method is similar to our case, while in the ADS method there are potentially 
even larger effects of order £d/A 2 and ?/£>/A 2 .) 

We note that while there are experimental advantages to using the resonant D decay 
final states K*K or pre, in general the full three-body Dalitz plot may be used to perform 
this analysis. (Higher multiplicity final states may also be used, although they will suffer 
from low reconstruction efficiencies.) Our method should also work well in some regions of 
the Dalitz plots of Cabibbo-allowed decay modes, such as Ks^n^ or KsK + K~ . However, 
the regions where and 1 — Tb are both of order unity are expected to be relatively small, 
and so the advantage of the high branching fractions of these dacays is not fully realized. 

In conclusion, we proposed a variant of the GW method to measure 7. It requires 
measurement of only Cabibbo-allowed and singly Cabibbo-suppressed D decays in color 
allowed B decays. Because it involves only large decay rates it may be carried out with 
current data sets. The branching fraction and reconstruction efficiency of the decay B ± — > 
K ± (KK*)e> is similar to that of B ± — > K ± (K + K~)d, which has already been observed jlll |. 
This might provide the first measurement of 7 that is free of hadronic uncertainties. 
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